Delivery of time-varying stimuli using ChR2 
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Abstract 



To understand sensory processing in neuronal populations, it is necessary to deliver stim- 
uli to the sensory organs of animals and record evoked population activity downstream. 
However, the pathways from sensory input to synaptic currents in cells that are several 
synapses removed from sensory organs are complex. Intrinsic noise and uncontrolled mod- 
ulatory input from other brain regions can interfere with the delivery of well-defined stim- 
uli. Here we investigate the ability of channelrhodopsins to deliver precise time-varying 
currents to neurons at any point along the sensory-motor pathway. To do this, we first 
deduce the amplitude response function of channelrhodopsin-2 (ChR2) using a three state 
Markov model of channel kinetics. With biophysically realistic parameters, this function 
supports a relatively broad signal passband and contains a resonance. We confirm the 
validity of our predicted amplitude response function using time-varying optical stimu- 
lation of excitatory neurons that express either wild type ChR2 or the ChR2(H134R) 
mutant. Together, our results indicate that ChR2-derived optogenetic tools are useful 
for delivering repeatable, time-varying currents to genetically-specified populations over 
a physiologically-relevant frequency band. 

Introduction 

The network response to time- varying input is at the core of cognitive and sensory process- 
ing. To understand how neuronal networks process time- varying input, precisely defined 
stimuli are delivered via a sensory organ, and evoked spiking activity is recorded down- 
stream. The power of this technique for deducing network encoding properties has been 
demonstrated in a number of sensory preparations, for instance, the retina and olfactory 
system [THH]. However, as stimuli delivered to sensory organs propagate to higher brain 
areas, intrinsic noise and modulatory input from secondary brain regions can interfere 
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with controlled input signals. For studies that seek to understand the function of neural 
circuits that are several synapses removed from sensory input, the direct introduction of 
time-varying currents to neural population may allow a more straightforward deduction 
of the circuit response properties. 

Optogenetic methods allow precise control of spike times using brief light pulses to 
excite light-gated cation channels, such as channelrhodopsin-2 (ChR2) [U(5]. Pulsed stim- 
ulation reliably dictates a spiking response that is tightly locked with each stimulus by 
briefly overriding neuronal dynamics. It stands in contrast to the highly variable, sub- 
threshold currents recorded from cortical neurons during natural sensory processing in 
vivo [6]. We hypothesized that using relatively low intensity, continuously modulated 
optical stimuli to excite ChR2 might allow conductance fluctuations that imitate natu- 
ral synaptic bombardment in vivo and leave the decision of when to spike to individual 
cells [TIE] - Surprisingly, while the response properties of microbial opsins to optical pulses 
have been studied extensively jl], little is known about their ability to relay fluctuating 
light signals. 

In order for ChR2 to be useful for delivering time-varying currents, it must allow 
(1) sufficient bandwidth to mimic synaptic communication and (2) repeatable current 
waveforms to be delivered across trials. Here, we address these requirements theoretically 
and experimentally. We find that wild-type ChR2 (wtChR2) [9] supports significant 
photocurrents from ~0-69 Hz and H134 mutant (ChR2(H134R) [10] from ~0-37 Hz. We 
show that evoked current waveforms are extremely repeatable across trials. Finally, we 
find that wild type ChR2 supports a strong resonance with a natural frequency around 
10 Hz. This resonance is present, but not significant in the H134R mutant. 
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Materials and Methods 



Derivation of ChR2's frequency response function 

To derive ChR2's light to photocurrent response function, we considered a three state 
Markov model of ChR2's state transition kinetics, which was introduced in [Tl] and is 
depicted in Fig. Ufa). The rate equations governing channel transitions are 

(1 ° (f ) -ecf>(t)C(t)-r d O(t) (i) 



dt 
dD(t) 



T d O(t) - r r D(t) (2) 



dt 

C(t) = l-0(t)-D(t), (3) 

where the state variables 0(t), D(t), and C(t) are the probabilities of a channel being open, 
desensitized, or closed, respectively. Yd and T r are the rates of channel desensitization 
and recovery, e is the quantum efficiency of ChR2 and 0(t) is the instantaneous photon 
flux for a single channel. <f>(t) can be modulated by changing the brightness of the optical 
stimulator. The evoked photocurrent across the cell membrane is proportional to the 
fraction of open channels. Therefore, under the assumption that a cell expresses a large 
number of channels, identifying ChR2's frequency response is equivalent to finding the 
frequency response of O in the continuum limit. We refer to ChR2's frequency response 
function as F ChR2 (u). 

i*ChR2(^) can be obtained by first considering a small sinusoidal light signal with a 
constant offset O , 

0(t) = 00 + 50 exp(jM), (4) 
where u = 2irf and / is the frequency of the sinusoid in Hz and j = y— 1. The first order 
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dynamics of the open and closed probabilities can then be described by a constant offset 
and small periodic component, 



0(t) = O + 5O exp(jcot) (5) 
D(t) = D + 5D exp(jwt). (6) 

In the linear regime, changes in the open state, 60, or the desensitized state, SD, are 
proportional to changes in optical input, e5<f). The proportionality factors for the open 
and desensitized states are the frequency response functions -FbhR2(w) and GchR2(w), 
respectively, 

50 = e50F ChR2 (u;) (7) 
SD = e500 ChR2 (u;). (8) 

Differentiating eqs. (J5J and (jSj) and inserting the result into eqs. and (J2J) leads to 

jtuSO exp(jcut) = [e0 o (l -0 - D )\ + [e50(l - O a - D Q ) (9) 

+ e0 o (5O - 5D) - T d 60\ exp(joot) + 0(2) 
ju5D exp(jwt) = [T d O - T r D ] + [T d 50 - T r 5D] exp(jut). (10) 

By dropping all but the first-order terms of eqs. (Q and (TTUT) (meaning those terms 
proportional to exp(jut)), and removing the common factor exp(jut), linear changes in 
the open and desensitized states due to changes in light power are given by 

joodO = e(l - O - D )6(f> + (e<j) - T d )50 - e<j)5D (11) 
judD = T d 50 - T r 5D, (12) 
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where (1 — O — D ) = C is the steady-state probability of the channel being closed. 
Performing the necessary substitutions to solve for 50 results in 



50 



C (juj + Y r ] 



-<jJ 2 + jco(Y r + e<p + Yd) + e<pY r + e<f)Y d + Y r Y d _ 



(13) 



Finally, referencing eq. fl7]), ChR2's frequency response function is calculated by dividing 
the left hand side of eq. f fT3"j) by e5(f), 

p ( U ) = *° = C,{ 3 u + Y r ) 

ChR2K } e5<P -u 2 +juj{Y r + e4) + Y d ) + e(t)Y r + e(t)Y d + Y r Y d ' K ' 



and the amplitude response is then given by 



"ChR2| 



V(-o; 2 + e0r r + e0r d + Y r Y d ) 2 + (oo(Y r + e0 + Y d )) 2 



(15) 



In the high frequency limit, eq. (fill) reduces to — 2 . c $") , , r , oc ■ /rr , c< i, r \ , i ■ This 
indicates that the frequency cutoff of i 7 bhR2(^) is determined by a low-pass filter with 
a time constant l/(r r + e0 + r<j). Thus increasing the baseline light intensity or the 
transition rates will result in a broader frequency range over which -FchR2(^) supports 
significant photocurrents. 

Experimental methods 

Culturing methods Our culturing methods are described in detail elsewhere [12] . All 
experiments were carried out in accordance with the U.S. Public Health Service's Policy 
on Humane Care and Use of Laboratory Animals and the Guide for the Care and Use 
of Laboratory Animals using a protocol approved by the Georgia Tech IACUC. Timed- 
pregnant female rats were anesthetized with inhaled isoflurane and killed by decapitation. 
Whole brains were excised from embryonic day 18 (E18) rats. Cortical tissue was digested 
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in a solution of 20 U-mU 1 papain (Sigma-aldrich). Following enzymatic digestion, cells 
were mechanically dissociated using 3 to 5 trituration passes through a plOOO pipette tip. 
The resulting cell suspension was filtered through with a 40 /im cell strainer and then 
centrifuged at 200-g to remove large and small debris, respectively. The cell pellet was 
diluted to 2500 cells-yuL -1 . Approximately 50,000 cells in a 20 /iL drop were plated at onto 
a ~2 mm diameter area on glass bottom petri dishes, resulting in ~2,500 cells-mm~ 2 on 
the culturing surface. 0.75 mL of the culturing medium was exchanged every three days, 
for each culture. Cultures dishes were sealed with a Teflon membrane [13] and stored in 
an incubator regulated to 35 °C, 5% CO2, 65% relative humidity. 

ChR2 expression system AAV2-CaMKlla::hChR2(H134R)-mCherry at 4-10 12 c.f.u.-mr 1 
was produced by the University of North Carolina at Chapel Hill Vector Core. AAV2- 
CaMKlla::wtChR2-mCherry at 4-10 12 c.f.u.-ml -1 was produced by the Kaplitt lab (Cor- 
nell University) using plasmid DNA for CaMKIIa::wtChR2-mCherry obtained from the 
K. Deisseroth (Standford University). At 1 to 5 days in vitro (DIV), viral aliquots were 
diluted to 1-10 12 c.f.u.-ml" 1 using glial-conditioned culturing medium. 1 fiL of diluted 
viral solution was added to 1 mL culturing medium for a final infection concentration of 
1-10 9 c.f.u.-ml -1 . Cultures were then incubated for 3 days before the culturing medium 
was exchanged. The fluorescent signal of the mCherry reporter protein was monitored 
in 3 sister cultures over the days post infection, and it increased monotonically before 
plateauing at ~3 weeks in vitro. All experiments were carried out on cultures that were 
3 to 4 weeks old. 

Optical stimulation A 10- watt light emitting diode (LED) was used for optical stim- 
ulation, with peak emission wavelength of 465 nm and ~20 nm full width at half maximum 
intensity (LZ4-00B200, LEDEngin, San Jose, CA). To deliver optical stimuli to cultured 
neurons, the LED was focused into the epi-illumination port of an E600FN upright mi- 
croscope (Nikon Corporation, Tokyo, Japan) and passed through a 40X objective lens. 
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The light power produced by LEDs is affected by their temperature. Additionally, the 
relationship between forward diode current and irradiance is a static non-linearity. To 
compensate for these factors and deliver distortion-free optical stimuli, we designed a pre- 
cision current source with integrated optical-feedback to drive our LED (Fig. |2^a)). This 
circuit measures the instantaneous optical power produced by the LED using an ampli- 
fied photodiode. It then adjusts the current sourced to the LED such that the optical 
power measurement matches a reference voltage supplied by a digital to analog converter 
(DAC; LIH 1600, HEKA Electronik, Lambrecht/Pfalz, Germany). The circuit can pre- 
cisely modulate the LED brightness over a bandwidth of 90 kHz (Fig. E^b-c)). A full 
design specification for the device is available online^. 

Intracellular recordings Whole-cell voltage-clamp recordings were conducted on 
neurons expressing the mCherry reporter protein. Recordings were performed in a con- 
tinuous perfusion of artificial cerebrospinal fluid (aSCF) bubbled with 95% 2 and 5% 
CO2 to maintain a pH of 7.4. The aSCF solution contained (in mM) 126 NaCl, 3 KC1, 2 
CaCl 2 , 1 NaH 2 P0 4 , 25 NaHC0 3 , 1.5 MgS0 4 and 25 D-glucose. The temperature of the 
extracellular medium was regulated to 35 °C using an inline heater (Warner Instruments, 
Hamden, CT). 1.5 mm outer diameter, 1.1 mm inner diameter borosilicate glass capilaries 
(Sutter Instruments, Novato, CA) were pulled into patch pipettes and filled with a solu- 
tion containing (in mM) 100 K-gluconate, 30 KC1, 3 ATP, 2 MgSO"4, 0.5 ethylene glycol 
tetraacetic acid and 10 HEPES adjusted to pH 7.4 using 0.1 M KOH. Filled pipettes had 
resistances of 4-8 M.Q. All recordings were performed using HEKA EPC8 amplifier and 
PatchMaster control software in whole-cell voltage clamp mode. Cells were held at -70 mV 
and membrane current measurements were amplified and low-pass filtered at 3 kHz be- 
fore being digitized at 20 kHz and streamed to disk. Access resistance and seal resistance 
were monitored between stimulation protocols. All experiments were performed in the 



1 https : //potterlab . gatech . edu/newman/wiki 
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presence of 40 /jM 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), 50 /iM (2R)-amino-5- 
phosphonovaleric acid (AP5), 20 [iM bicuculline in order to prevent synaptic transmission. 
Whole-cell recordings were analyzed offline in MATLAB (The MathWorks, Natick, MA). 



Experimental verification of frequency response functions 

To estimate the frequency response of ChR2, we delivered optical stimuli, s(t), consist- 
ing of T=10 second realizations of a Gaussian (Ornstein-Uhlenbeck) noise process while 
recording evoked intracellular currents, h{t), within a single cell, i. s(t) was generated 
using 



where s(ti) — 0, [x s — 0.4, and a s — 0.08 mW-mm -2 are the initial condition, mean, and 
standard deviation of the process respectively, r s = 50 ms is the correlation time s(t), 
dt = 40 [is is the DAC update period, and £(t n ) is a random variable drawn from the 
standard normal distribution. Each cell was exposed to the a single, repeated realization 
of s(t) for k = 10 trials. The first second of each trial was ignored to remove the non- 
stationary effects of the stimulator turning on. The recorded intracellular currents were 
averaged across trials, 



k=i 

to remove trial-to-trial noise. We then used reverse correlation techniques to estimate the 
empirical frequency response function for each cell, 



s(Wi) = A*s + s(t n )exp(-dt/r s ) + a 3 ^(t n )y/l - exp(-2rft/r s ), 



(16) 




(17) 



^ChR2,i(w) 



S s(h) 



(18) 
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where S ss is the power spectrum of s(t) and S' s /j i \ is the cross spectrum of (Jj) and s(t). S ss 
and Sgfa) are defined as the Fourier transforms of the corresponding correlation function, 



Finally, we averaged FchR2,i( u ) across ce hs to obtain the empirical frequency response 
for each construct, -FchR2(w). To improve our estimate of the power spectra, we followed 
the procedure introduced in [TJ] and used a frequency dependent window, equivalent to a 
Gaussian bandpass filter with standard deviation of o = 2n/u in the frequency domain. 
Spectra were evaluated at discrete increments, u n = 27rl0™, n = 0.1, 0.2, 3 . 



In this study, we sought a general description of ChR2's dynamics that captured the ability 
of both wild type ChR2 (wtChR2) and its engineered variants to convey time-varying 
stimuli. To do this, we determined the frequency response function of a population of 
channels expressed by a single cell, i 7 chR2(^) ; using a Markov model of ChR2's channel 
kinetics [TTJ (see Methods for details). The model is illustrated in Fig. [T](a) and F ChR2 {^) 
is given by eq. ffT4l) . The current that is passed by ChR2 across the cell membrane is 
proportional to the number of channels that occupy the open state. Therefore, -Fcm^O^) 
can be thought of as a frequency- and phase-dependent transition rate from the channels' 
closed to open state in response to a time-varying stimulus. Since individual channels 
switch between states discretely, -FbhR2(w) describes the transfer of optical waveforms to 
intracellular current under the assumption that a large number of channels are present in 
the cell's membrane. 




(19) 



(20) 



Results 
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The amplitude response function, \F wt c hR2 (u)\, eq. (TIB]), indicates the frequency depen- 
dent gain of the channel population to incoming light signals. The predicted amplitude 
response function for wtChR2, \F wt chR2{u)\, and the H134R mutant, |i 7 bhR2(Hi34R)(w)|, 
are shown in Fig. [H(b) for different mean illumination intensities. |-F w tchR2(^)| has high 
frequency cutoff (full width at half maximum) of 69 Hz. It should be noted that this cutoff 
value is defined relative to wtChR2's peak conductance, and not in terms of absolute pho- 
tocurrents. For this reason, it is possible to use wtChR2 to deliver physiologically signifi- 
cant photocurrents at stimulation frequencies exceeding 69 Hz. The shape of |-F w tchR2(w)| 
indicates that wtChR2 exhibits a significant resonance with a natural frequency close to 
10 Hz. This feature explains the large peak to steady-state (DC) current ratio exhibited 
by wtChR2 [4j. In agreement with previous characterizations of the channel, the H134R 
variant is significantly slower than wtChR2 and |i 7 bhR2(Hi34R)(w)| has a cutoff frequency 
at 37 Hz. While ChR2(H134R) supports a resonance, it is significantly reduced compared 
to wtChR2. 

To verify -F ChR2 (w) experimentally, cultured cortical cells expressing either wtChR2 
or ChR2(H134R) were stimulated with 465 nm at peak intensity, spatially uniform light 
while somatic photocurrents were recorded using whole cell patch clamp (see Methods for 
details). Fig. [2]^a-c) details the stimulation system used in our study. Optical stimuli con- 
sisted of a 10 second realization of a Gaussian noise process with a time constant of 50 ms 
and a mean±STD irradiance of 0.4±0.08mWmm -2 . We choose stimuli with these param- 
eters because they evoked currents with similar amplitude and temporal characteristics 
to those obtained from in vivo recordings of sensory cortical neurons [6l[T5l[T6] . 

We estimated the empirical frequency response function, -F ChR2 (w), for cells expressing 
wtChR2 (n = 9 cells) or ChR2(H134R) (n = 4 cells) using reverse correlation techniques 
(eq. (1T81) ). Fig. [T]^c) compares the empirical amplitude response, |-fchR2(w)|, of wtChR2 
and ChR2(H134R) and their theoretical counterparts. We observed good agreement be- 
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tween the predicted response function and |FchR2(^)|, although some differences exist. 
For instance, both |i 7 chR2(Hi34R)(^)| and \F wtC hR2(u)\ have a small downward deviation 
from the predicted response at ~5 Hz, which is more pronounced for wtChR2. Addition- 
ally, the predicted frequency response tends to slightly overestimate the measured gain 
at frequencies above 100 Hz. Despite these differences, both theoretical and empirical 
results indicate that wtChR2 and ChR2(H134R) are useful tools for transmitting fluctu- 
ating current stimuli to populations of cells in a physiologically relevant frequency range. 
Additionally, because the model provides an extensible description of channel dynamics, 
it serves as a useful tool for predicting the bandwidth and resonance of new channels 
based on measurable physiological parameters. 

Next, we measured the reliability of photocurrents across trials. Fig. [21(d) shows a 
3-second portion of a stimulus waveform and the corresponding photocurrents for a single 
cell expressing wtChR2, across trials. As expected, these photocurrents appear to be 
smoothed versions of the stimulus signal due to the low-pass effect of wtChR2's amplitude 
response function. Photocurrents were remarkably stable across trials. The trial-to-trial 
repeatability of evoked photocurrent waveforms is captured in Fig. EJ^f), which displays 
the standard deviation of the current waveform recorded during the first presentation 
of the stimulus compared to the 10th presentation, for each cell. Because these points 
fall near the identity line, it can be inferred that there is no change in the efficacy of 
photostimulation across trials. 

We then examined the distribution of evoked current amplitudes across cells ( Fig.[2](f)). 
wtChR2 delivered current waveforms with a mean waveform standard deviation of 26.7 pA 
per cell. ChR2(H134R) delivered only slightly larger current fluctuations than wtChR2, 
with a mean waveform standard deviation of 32.0 pA per cell. The comparable amplitudes 
of evoked currents between wtChR2 and ChR2(H134R) is due to wtChR2's strong reso- 
nance, which makes it more sensitive to fluctuating, compared to DC, optical waveforms. 
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Finally, to determine the reliability of evoked currents across cells, we calculated the 
normalized cross- correlation function, c^, between the light power density, s(t) and pho- 
tocurrents Ii(t) at each cell, i, and across cells, ci i} i r Fig. [I] (g) and (h) show c s ^ and 
cj u i. for wtChR2. The median peak value of cij j is 0.96, indicating strong correlations 
between evoked currents in different cells. The median peak value of c s i. is 0.92, indicating 
strong correlations between evoked currents and the stimulus waveform. Additionally, the 
similarity in shape between c s j and the autocorrelation of the stimulation process, c ss , 
indicates that precise temporal features of the stimulus were accurately converted into a 
photocurrent, as predicted by the passband of the frequency response functions. 

Discussion 

Optogenetic methods offer genetic specificity, elimination of electrical recording artifacts, 
and an expanding toolset with increasingly specialized functionality [4]. Because of these 
advantages, optogenetic methods are often used for direct interaction with neuronal sub- 
populations in order to understand their function [TTt[T8]. Typically, pulsed optical stimuli 
are used to probe neural circuits. However, for studies that seek to understand informa- 
tion transmission in neural circuits, continuously modulated photocurrents better mimic 
natural synaptic bombardment than do pulsed stimuli. 

In this study, we demonstrated the ability of ChR2 to deliver continuously modulated 
photocurrents to neurons. We used a three state Markov model [H] to derived an an- 
alytical frequency response function for ChR2 variants (eq. ( fl4l) ). We confirmed this 
model experimentally indicating that eq. ( fl4|) is sufficient to capture dynamical prop- 
erties of ChR2 in neurons. Additionally, we found that the passband of wtChR2 and 
ChR2(H134R) is broad enough to support photocurrents that mimic natural synaptic 
communication. The reliability of somatic photocurrents is very high, with a correla- 
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tion strength upward of 0.9 across trials and cells. ChR2's amplitude response function, 
eq. (115p . indicates that the sum of channel recovery and desensitization transition rates 
determines the frequency cut-off. Therefore, opsins with faster transition rates [Sj will 
likely allow an even broader dynamic range. Interestingly, faster channels and increased 
bandwidth may eventually offer an artificial, optical neural communication channel that 
actually exceeds the bandwidth of natural sensory organs. This would have tremendous 
implications for how neural computation and processing are studied and for the advance- 
ment of brain-machine interfaces. Finally, we showed that ChR2's frequency response 
function supports a resonance. The degree of resonance is dependent on the values of 
free model parameters, which will change for different ChR2 variants. This finding is 
especially important for studies that use ChR2 to examine the frequency-dependence of 
neural circuitry [17], since it is important to separate the intrinsic dynamics of ChR2 from 
those that belong to the network under study. 

ChR2 was derived from microbes that use it for optical sensation in natural environ- 
ments. It is therefore not surprising that the channel is excellent at conveying time- varying 
input signals. Using wtChR2 and its numerous variants as a means for delivering repeat- 
able, time-varying stimuli to genetically defined populations of cells will be a powerful 
method for probing the dynamics of neural circuits. 
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Figures 
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Figure 1: ChR2's amplitude response function, (a) Illustration of the three-state Markov channel 
model described by eqs. (JJ) - The transition rates between open, O, desensitized, D, and closed, 
C, states are ecj)(t), T r , and IV (b) Amplitude response functions for the model, eq. (fl~5]). is shown for 
different mean illumination intensities (0.15 through 0.6 mW-mm~ 2 ). (c) Predicted amplitude response 
of wtChR2 (solid lines) and ChR2(H134R) (dashed lines) compared to experimentally measured response 
functions for a mean illumination intensity of 0.4 mW-mm -2 . Error bars are ±1 STD (n = 9 neurons for 
wtChR2 and n = 4 neurons for ChR2(H134R)). 
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Figure 2: Reliability of time-varying photostimulation. (a) Simplified schematic of LED driver 
in optical feedback mode, (b) One millisecond LED pulse, Vpd (black), versus the reference voltage, 
Vref (gray). Inset shows the step onset with corresponding scale bars. The current sourced to the 
LED is shown in the bottom plot, (c) A DAC-controlled Gaussian stimulus (gray) signal and recorded 
light waveform (black). Inset shows zoomed portion of the sequence. An amplitude histogram of the 
sequence, with best-fit Gaussian, is shown to the right, (d) Intracellular currents, from a single cell, 
during Gaussian stimuli. The top trace is a portion of a 10 second Gaussian stimulus sequence. The 
bottom three traces show the intracellular currents recorded during different presentations of the same 
stimulus waveform to a single cell, (e) The same stimulus waveform used in (d), and the corresponding 
evoked responses from different cells, (f) The standard deviation of the photocurrent induced on the 
first trial of stimulation versus the last trial. The dotted line is identity. The filled dot corresponds to 
the cell in (d). (g) Cross/auto-correlation functions of photocurrents between neurons, Ii(t). The inset 
histogram shows the peak correlation coefficients, (h) Cross-correlation function between the stimulation 
process s(t) and recorded photocurrents. The grey line is the autocorrelation function of the stimulation 
process. The inset histogram shows the peak correlation coefficients. Unless otherwise noted, data in this 
figure were obtained from cells expressing wtChR2. 
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